Surface Partition of Large Fragments 
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The surface partition of large fragments is derived analytically within a simple statistical model 
by the Laplace-Fourier transformation method. In the limit of small amplitude deformations, a 
suggested Hills and Dales Model reproduces the leading term of the famous Fisher result for the 
surface entropy with an accuracy of a few percent. The surface partition of finite fragments is 
discussed as well. 



I. INTRODUCTION 

During last forty years the Fisher droplet model 
(FDM) y , on one hand, was extensively used to analyze 
the condensation of gaseous phase (droplets of all sizes) 
into hquid. The gaseous phases are ranging from mixture 
of nuclear fragments i2j to the various clusters on 2- and 
3-dimensional Ising lattices Q- On the other hand, the 
FDM inspired the formulation of a more sophisticated 
model, the statistical multifragmentation model (SMM) 
U , which describes not only the gaseous phase of nuclear 
fragments, but the liquid phase (nuclear matter) as well 
on the same footing ^ l^j Q • 

Applying the FDM to the nuclear fragment of A- 
nucleons, one can cast its free-energy Fa as follows 

Fa^-W A + a{T) A^/^ +tT In A (1) 

Here W is the bulk binding energy per nucleon, cr(T) is 
the temperature dependent surface tension which about 
the critical temperature Tc is parameterized in the fol- 
lowing form: a{T) ^ ao[l - T/Tc], with CTo ~ 18 MeV 
and Tc = 18 MeV (ct = at T > T^). The last con- 
tribution in Eq. ^ generates the famous Fisher's term 
with dimensionless parameter t. From the study of the 
combinatorics of lattice gas clusters in two dimensions, 
Michael Fisher had postulated Eq. (^ and this kind 
of temperature dependence of the surface tension be- 
cause the latter naturally leads to the existence of crit- 
ical temperature Tc- This is, of course, not a unique 
parametrization of the surface tension. The SMM, for 
instance, successfully employs another one 4] a{T) = 
cro[(T2 - T^)/{T^ + T2)]5/4. Therefore, it is necessary 
to study a few simple, but quite fundamental questions, 
"What is the origin of the Fisher parametrization for the 
temperature dependent surface tension in three dimen- 
sions? Does any physical motivation favor the Fisher 
parametrization?" This work is devoted to these ques- 
tions of fundamental importance. 



II. HILLS AND DALES MODEL 

To answer these main questions we will consider the 
statistical model of surface deformations. We will, how- 
ever, impose a necessary constraint that the deforma- 
tions should conserve the total volume of the fragment 



of A-nucleons. As we will see, the most interesting result 
corresponds to the deformations of vanishing amplitude. 
Therefore, it is clear that the shape of deformation can- 
not be important and we can choose one which is regular 
to simplify our presentation. For this reason we shall con- 
sider cylindrical deformations of positive height hk > 
(hills) and negative height —hk (dales), with k nucleons 
at the base. For simplicity it is assumed that the top 
(bottom) of the hill (dale) has the same shape as the 
surface of original fragment of A nucleons. Our main 
assumptions are as follows: (i) the statistical weight of 
deformations exp (— tTo|A5'fc|/T/si) is given by the Boltz- 
mann factor due to the change of the surface |AS'fc| in 
units of the surface per one nucleon si; (ii) the hill's 
heights hk < Hk {Hk is the maximal height of the hill 
with k nucleons at the base) have the same probability 
dhk/Hk besides the statistical one; (iii) the same two as- 
sumptions are valid for the dales as well. As we will see, 
these assumptions are not too restrictive, but allow us to 
simplify the analysis. 

Under adopted assumptions it is possible to find the 
one-particle statistical partition of the deformation of the 
k nucleons base as a convolution of two probabilities: 
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where upper (lower) sign corresponds to hills (dales). 
Here Pk is the perimeter of the cylinder base. Now we 
have to find a geometrical partition (degeneracy factor) 
or the number of ways to place the center of a given de- 
formation on the surface of A-nucleon fragment while it 
is occupied by the set of {ri; = 0, 1, 2, ...} deformations 
of the / nucleons base. Our next assumption is that the 
desired geometrical partition can be given in the excluded 
volume approximation 



Sa- Y. ^i'^k + "fe )*i 
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Si 



(3) 



where sik is the area occupied by the deformation of k 
nucleon base (k — 1,2, ...), Sa is the full surface of the 
fragment, and KmaxiSA) is the A-dependent size of the 
maximal allowed base on the fragment. It is clearly seen 



now that the denominator in the right hand side (r.h.s.) 
of © corresponds to the available surface to place the 
center of each of {n^} deformations that exist on the 
fragment surface. It is necessary to impose the condition 
G > which ensures that the deformations do not over- 
lap on the available surface of the fragment. Eq. JSJ is 
the Van der Waals excluded volume approximation usu- 
ally used in statistical mechanics at not too high particle 
densities [3,11, 111 and it can be derived for the objects of 
different sizes in a spirit of a method proposed in [9( . 

According to Eq. ^ the statistical partition for the 
hill with a fc-nucleon base matches that of the dale, i.e. 
■^k ~ ^k ■ Therefore, the grand canonical surface parti- 
tion (GCSP) 
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corresponds to the conserved (on average) volume of the 
fragment because the probabilities of hill and dale of the 
same base are identical. The presence of 0(sit/)-function 
in I0J ensures that only configurations with positive value 
of the free surface of fragment are taken into account. 
However, this makes the calculation of the GCSP very 
difficult. Because of the explicit Sa dependence of the 
maximal base of deformations via Kmax{SA) the stan- 
dard trick to deal with the excluded volume partitions, 
the usual Laplace transform method _5^, _6, _8| in Sa^ can- 
not overcome this difficulty. However, the GCSP Q can 
be solved with the help of the recently developed Laplace- 
Fourier technique ^]. The latter employs the identity 
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which is based on the Fourier representation of the Dirac 
J-function. The representation (O allows us to decouple 
the additional ^A-dependence in Kmax{SA) and reduce it 
to the exponential one, which can already be integrated 
by the Laplace transform pil | 
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After changing the integration variable Sa — > 5' = Sa ~ 

^ k (n^ + rij. ) Si, the constraint of 6-function has 
fe=i 



disappeared. Then all n^ were summed independently 
leading to the exponential function. Now the integration 
over S' in ^ can be straightforwardly done resulting in 
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where the function J-'{£,, A) is defined as follows 
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As usual, in order to find the GCSP by the inverse 
Laplace transformation, it is necessary to study the struc- 
ture of singularities of the partition |7J). Since the HDM 
requires the fragment volume conservation, hereafter we 
will call ((TJ as an isochoric partition or isochoric ensem- 
ble. 



III. ISOCHORIC ENSEMBLE SINGULARITIES 

For finite fragment surface the structure of singularities 
of the isochoric partition ((JJ) can be complicated. To see 
this let us first make the inverse Laplace transform: 
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where the contour integral in A is reduced to the sum over 
the residues of all singular points X = Xn + irj with n = 
0, 1, 2, .., since this contour in the complex A-plane obeys 
the inequality x > max(i?e{A„}). Now all integrations 
in ^ can be done, and the GCSP acquires the form 
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i.e. the double integral in © simply reduces to the sub- 
stitution ^ —f Sa in the sum over singularities. This 
remarkable answer is a partial example of the general 
theorem on the Laplace-Fourier transformation proper- 
ties proved in \W\ . 

The simple poles in are defined by the condition 
Xn — J-{Sa, Xn) and the latter can be cast as a system 



of two coupled transcendental equations 

Rn= J2 H + ^k] e-^-«"cos(/„/c), (11) 

In = - E H + ^k] e-^-^"sm(/„fc), (12) 
fc=i 

for dimensionless variables i?„ = si_Re(A„) and /„ = 
si/?Ti(A„). So far Eqs. (|11() and H12|l are rather gen- 
eral and can be used for particular models which specify 
the height of hills and depth of dales. But we can give 
an absolute supremum for the real root {Rq]Iq — 0) of 
these equations. For this purpose it is sufficient to con- 
sider the limit Kmax{SA) ~^ oo, because for /„ = /q = 
the right hand side (r.h.s.) of (|ll|l is a monotonously 
increasing function of Kmax{SA)- Since z^ = z'^ are the 
monotonously decreasing functions of Hk, the maximal 
value of the r.h.s. of I|1H) corresponds to the limit of in- 
finitesimally small amplitudes of deformations, Hk — > 0. 
Under these conditions Eq. (|12|l for /„ = /q = becomes 
an identity and Eq. (|11|1 acquires the form 
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and we have Rq ~ siXq w 1.06009. Since for /„ j^ 
defined by H12() the inequality cos(/„fc) < 1 cannot be- 
come the equality for all values of k simultaneously, then 
it follows that the real root of pi|) obeys the inequality 
Rq > Rn>o- The last result means that in the limit of 
infinite fragment, Sa —* oo, the GCSP is represented by 
the farthest right singularity among all simple poles {A„} 



Z{Sa) 



Sa- 



e si 
l + # 



RqSa 

0.6402 e «i 



(14) 



There are two remarkable facts regarding (|14|l : first, this 
result is model independent because in the limit of van- 
ishing amplitude of deformations all specific parameters 
of the model have dropped out; second, in evaluating 
(|14|l we did not specify the shape of the fragment un- 
der consideration, but only implicitly required that the 
fragment surface together with deformations is a regular 
surface without self-intersections. Therefore, for vanish- 
ing amplitude of deformations the latter means that Eq. 
p4|l should be valid for any self-non-intersecting surfaces. 
For spherical fragments the r.h.s. of H14|l becomes fa- 
miliar, 0.6402 e-*^-^"^^^^ , which, combined with the 

Boltzmann factor of the surface energy e ° ' , 
generates the following temperature dependent surface 
tension of the large fragment 
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FIG. 1: The first quadrant of the complex plane siA,i = 
Rn+iln of the roots of the system of (llfcil and 11711 . The circles 
and squares represent the two branches I^ and I^ of the 
roots, respectively. The curve is defined by the approximation 
given by 11911 and 1201 1 (see text for more details). 



which means that the actual critical temperature of 
the three dimensional Fisher model should be Tc ~ 
(To /1. 06009, i.e. 6 % smaller than Fisher originally sup- 
posed. Note please that this equation for the critical 
temperature remains valid for the temperature depen- 
dent (To as well. It is also surprising that the degeneracy 
factor 0.6402 in (|14|l is only 12.5 % larger than the cor- 
responding factor of the self-avoiding polygons on the 
two-dimensional square lattice |ll| . 

For the large, but finite fragments it is necessary to 
take into account not only the farthest right singularity 
Ao = Ro/si in l(Tn|l . but all other roots with positive real 
part i?,i>o > 0. In this case for each i?„>o there are two 
roots ±/„ of (|12|l because the GCSP is real by defini- 
tion. The roots of Eqs. pi|l and H12I) with largest real 
part are very insensitive to the large values of KmaxiSA), 
therefore, it is sufficiently good to keep Kmax{SA) -^ oo. 
Then for limit of vanishing amplitude of deformations 
Eqs. Hll|) and (|12|l can be, respectively, rewritten as 
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After some algebra the system of (|16f) and (|17|) can be 
identically reduced to a single equation for Rn 
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and the quadrature /„ 
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— i?^j. The analy- 



sis shows that besides the opposite signs there are two 
branches of solutions, /+ and /~, for the same n > 1 
value: 
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The exact solutions {Rn',In) for n > 1 which have the 
largest real part are shown in Fig. 1. together with 
the curve parametrized by functions I^ and Rx taken 
from Eqs. lfTO|l and (|2DJ), respectively. From Eq. l|^ 
and Fig. 1. it is clearly seen that the largest real part 
i?i w 0.0582 is about 18 times smaller than Rq, and, 
therefore, already for the fragment of a few nucleons the 
correction to the leading term H14() is exponentially small. 
Using the approximations 119|) and 120|) . for n > 2 one 
can estimate 
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the upper limit of the (i?„ ; /j^ ) root contribution into the 
GCSP IjlOl) . This result shows that the total contribution 
of all complex poles in IjlOl) is negligibly small compared 
to the leading term (|14|l for a fragment of a few nucleons 
already. The latter, however, requires a more careful 
accounting for the volume conservation of a fragment. 



IV. CONCLUDING REMARKS 

The developed model allows us to give the upper limit 
for the surface entropy because it corresponds to the van- 
ishing amplitude of deformations and, therefore, the spe- 
cific features of the model were irrelevant for our analysis. 
To find the next order correction to the surface entropy 
one has already to consider the underlying model for de- 
formations. We, however, would like to show how the 
power law may arise within the HDM. For this purpose 



let us consider the left equality in (|13|l which is valid for 
small heights of deformations. It can be shown that the 
following ansatz (Sa >> si) for the deformation energy 
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of k nuclcon base, indeed, generates the Fisher power law 
A^'^ for the GCSP H14|l of an A-nucleon fragment. From 

(|22|l it is clearly seen that the term — fc JJ^ In ( 4^ j is the 

entropy which gives an a priori uncertainty to measure 
the position of k nucleons of area si on the surface of the 
fragment. A comparison of (|22|l with any fci?„ > in the 
left equality H13(l shows that in the limit 5*^ >> si the 
ansatz (I22|l corresponds to a negligible correction com- 

pared to the exponentials e =i . Therefore, the Fisher 
power law is too delicate for the model developed to 
study the surface partition of large fragments. 

In conclusion, we developed a simple statistical model 
which allowed us to derive analytically the general ex- 
pression (|10|l for the GCSP. This result is achieved by 
the Laplace-Fourier transformation technique to the iso- 
choric ensemble. We named this ensemble because the 
HDM obeys the volume conservation of a fragment un- 
der consideration. The volume conservation is accounted 
for by the equal statistical probabilities for the hills and 
dales of the same base. The limit of vanishing deforma- 
tions allowed us to find a supremum for the surface en- 
tropy of the large fragments which, remarkably, exceeds 
the Fisher original assumption by only 6 percent. The 
analytical analysis of the corrections to the GCSP H14|) 
originating from the complex roots of Eqs. (|ll|l and H12|l 
showed that these corrections are negligible already for 
the fragment of a few nucleons. The HDM allows one to 
study the statistical mechanics of volume deformations 
of finite or even small fragments, but this task requires 
further refinements of the model. 
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